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Charged current inclusive neutrino-nucleus cross sections are evaluated using the superscaling 
model for quasielastic scattering and its extension to the pion production region. The contribution 
of two-particle-two-hole vector meson-exchange current excitations is also considered within a fully 
relativistic model tested against electron scattering data. The results are compared with the inclusive 
neutrino-nucleus data from the T2K and SciBooNE experiments. For experiments where (E v ) ~ 
0.8 GeV, the three mechanisms considered in this work provide good agreement with the data. 
However, when the neutrino energy is larger, effects from beyond the A also appear to be playing 
a role. The results show that processes induced by two-body currents play a minor role at the 
kinematics considered. 
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I. INTRODUCTION 

New measurements of inclusive charged current (CC) 
neutrino-nucleus scattering cross sections, where only the 
outgoing lepton is detected, have been recently performed 
by the T2K 00, SciBooNE i and ArgoNeuT 0,1 col¬ 
laborations. For neutrino energies around 1 GeV (T2K 
and SciBooNE) the main contributions to the cross sec¬ 
tions are associated with quasiclastic (QE) scattering and 
one pion (l7r) production. These, along with the two- 
particle-two-hole (2p2h) meson-exchange current (MEC) 
contributions, are the only processes we shall consider in 
this paper, while at higher neutrino energies (ArgoNeuT) 
multiple pion and kaon production, excitation of res¬ 
onances other than the A and deep inelastic channels 
should also be considered. That said, we do attempt to 
provide some insights into how important these last ef¬ 
fects may become as the neutrino energy increases. 

The QE muon neutrino and antineutrino cross sections 
measured by the MiniBooNE experiment 0 7}, where QE 
events are characterized by the absence of pions in the 
final state, have triggered a lot of theoretical work try¬ 
ing to explain the unexpectedly large results, in appar¬ 
ent tension with the higher-energy data from the NO¬ 
MAD experiment [§}. Several calculations [Ml have 
demonstrated that 2p2h excitations induced by two-body 
meson-exchange currents (MEC) play a significant role in 
the interpretation of the QE MiniBooNE data and in the 
neutrino energy reconstruction, which is therefore model 


dependent. At a quantitative level, however, these calcu¬ 
lations, relying on different models and approximations 
(see [131 ] for a brief review of the various approaches), 
give quite different results. Furthermore, a model using 
Relativistic Green’s Functions (RGF), which does not ex¬ 
plicitly contain two-body currents but to some extent in¬ 
cludes inelastic channels through a complex optical po¬ 
tential, has been shown to be able to explain the QE 
MiniBooNE data fl4| : however, these results depend sig¬ 
nificantly on the particular choice of the relativistic op¬ 
tical potential. Comparison with inclusive data, where 
many-particle (in particular two-nucleon) emission chan¬ 
nels unambiguously contribute, can shed light on the role 
of MEC in neutrino and antineutrino scattering in differ¬ 
ent kinematical regions. 

In this paper we evaluate the CC neutrino inclusive 
cross sections within the Superscaling approach (SuSA), 
introduced in [15] to describe neutrino-nucleus scatter¬ 
ing by using electron scattering data instead of relying 
on specific nuclear models. This approach allows one to 
describe the QE and A resonance regions in a unified 
framework and can be applied to high energies due to its 
relativistic nature. In the QE region, the SuSA model has 
recently been improved in [16j to incorporate effects aris¬ 
ing in the Relativistic Mean Field (RMF) model in the 
longitudinal and transverse nuclear responses, as well as 
in the isovector and isoscalar channels. Since MEC are 
known to violate superscaling, their contribution must 
be added to the superscaling result: this will be accom- 
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plished by using a parametrization pa of the relativistic 
calculation of [18l |. The parametrization is necessary in 
order to reduce the computation time, since the exact 
calculation involves 7-dimensional integrals. 

The paper is organized as follows: in Sect. Ql] we briefly 
describe the Superscaling approach to the QE and A- 
resonance regions, including some recent improvements 
of the model. In Sect. m we test our approach with 
inclusive electron scattering data on 12 C and we present 
the comparison of the calculation for muon-neutrino and 
electron-neutrino scattering with the data from the T2K 
and SciBooNE experiments; in the latter case we also 
present results for antineutrinos. In Sect. m we draw 
our conclusions. 


II. THE SUPERSCALING MODEL 


The SuperScaling model, based on the superscaling 
properties of inclusive electron scattering EH, |20|, has 
been extensively used ( 2 TI [25^ to predict neutrino and an¬ 
tineutrino cross sections for complex nuclei. The detailed 
description of the model can be found, e.g., in 00 . 
Here we simply recall its main features. In the quasielas¬ 
tic peak (QEP) region the basic ingredient of the model 
is a phenomenological superscaling function 
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extracted from the world electromagnetic (e,e') data by 
dividing the longitudinal response R^ E times the Fermi 
momentum kp by the single-nucleon elementary function 
G® E . The data show that f® E is to a large extent a 
function of only one variable, the scaling variable ipQp, 
and is independent of the momentum transfer q (scaling 
of first kind) and of the nucleus, represented by the Fermi 
momentum ftp (scaling of second kind). 

The function f® E embeds most of the nuclear effects, 
both in the initial and in the final state, and can therefore 
be used to predict the weak charged current quasielastic 
(CCQE) (vi,l) cross section. In its original version the 
SuSA model assumes that the superscaling function f® E 
is the same in the longitudinal (L) and transverse (T) 
channels, a property referred to as scaling of zeroth-kind. 

The main merit of the SuSA model is the reasonable 
agreement, required by construction, with electron scat¬ 
tering data over a very wide range of kinematics and for 
a large variety of nuclei. Such an agreement is a cru¬ 
cial test for any nuclear model to be applied to neutrino 
reactions. Although phenomenological, the model has 
firm microscopic foundations in the RMF model, which 
is able to reproduce both the height and the asymmet¬ 
ric shape of the experimental superscaling function [ 2 (| . 
Furthermore, and importantly, the RMF model predicts 


a transverse superscaling function 

jdQE 

f? E = k F % ( 2 ) 
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which is higher than the longitudinal one, a result sup¬ 
ported by the separated L/T data analysis EH, H3 and 
strictly linked to the relativistic nature of the model [28| . 
This result has recently been used to improve the ingre¬ 
dients of the SuSA model by constructing a new version 
(SuSAv2) where f® E > EH- Moreover in SuSAv2 
the effects of Pauli blocking, initially neglected, have been 
implemented. In the results we present in the next sec¬ 
tion the updated version SuSAv2 of the model will be 
used. 

The scaling approach can then be inverted and pre¬ 
dictions can be made for CCQE neutrino and antineu¬ 
trino reactions by replacing the elementary electromag¬ 
netic vertex, 7 *NN , with the weak one, WNN. 

In the QE region the superscaling predictions have 
been succesfully compared with the recent MINERVA 
data [29|, H 3 | , that have been shown in EMI to be well 
reproduced without need of invoking large 2 p 2 h contri¬ 
butions. Good agreement is also obtained with the high- 
energy NOMAD data [§[. On the contrary, the Mini- 
B 00 NE QE data are underpredicted by the model. The 
inclusion of 2p2h MEC excitations in the vector channel, 
evaluated using the model of EH , gives results which are 
closer to the experimental points, but are not enough to 
explain the data Enm, unlike the analysis and results 

in mm- 

It should be mentioned that the model developed 
in EH f° r electron scattering only contains the vector 
part of the two-body current. Assuming the transverse 
vector 2p2h MEC scaling function, fjfyy, to be equal 
to the axial-axial (f-pj fjf) and vector-axial (fp va) ones 
— as considered in |9| — a final result in agreement with 
CCQE MiniBooNE data is found. However, such a result 
cannot be fully justified until a proper 2p2h MEC cal¬ 
culation for the axial-axial and vector-axial responses is 
completed. The full calculation, including the axial two- 
body current, is in progress [34[ and once it is available 
it will allow us to test the quality of this approximation. 

The superscaling approach has been extended from the 
QE domain into the region where the A-excitation dom¬ 
inates. In EH, mi it has been shown that the residual 
strength in the resonance region, obtained by subtract¬ 
ing the QE contribution from the total cross section, can 
be accounted for by introducing a new scaling function 
f A dominated by the N —¥ A and employing a new scal¬ 
ing variable, , which is suited to the resonance region. 
In this paper we revisit this approach by using the im¬ 
proved QE superscaling model, SuSAv2, and an updated 
parametrization of the 2p2h MEC response EH- This 
procedure yields a good representation of the electromag¬ 
netic response in both the QE and A regions, as we shall 
illustrate in the next section. 

Two different parameterizations have been considered 
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to deal with the less-known axial form factors C 3 4 5 ap¬ 
pearing in the elementary W + N —» A + transition cur¬ 
rent. One is taken from [36( where the deuteron was stud¬ 
ied and the other was introduced in [13] • The comparison 
of results obtained with the two parameterizations can 
be viewed as a measure of the degree of uncertainty that 
can be expected from the choice of the single-nucleon 
response for this reaction. In the present analysis, our 
results show a negligible dependence upon the choice of 
parametrization. Hence all results presented in this pa¬ 
per correspond to the model of [§(| • 

The superscaling predictions for the Delta region have 
been compared in [38} with the MiniBooNE data in the 
case of 7 r + production for the v^-CB .2 CC charged pion 
cross section [13 ■ The results obtained for the flux- 
averaged double- and single-differential cross sections as 
functions of the muon kinetic energy and scattering an¬ 
gle were found to be in good agreement with the data, 
whereas for the totally integrated unfolded cross section 
a somewhat different dependence on the neutrino energy 
was obtained from the one displayed by the data. It is 
also important to stress that the present scaling approach 
is expected to be valid only for those kinematical situa¬ 
tions where the A-resonance excitation is the dominant 
inelastic process. At higher energies heavier resonances 
can be excited and the deep inelastic scattering domain 
can be reached. In this case the phenomenological exten¬ 
sion of the model developed in J40j> based on direct fits 
to highly inelastic e — N scattering data, is more suitable 
to describe the inclusive cross section. 

Finally, as first noticed in (l5) and studied in more 
depth in [35}, sizeable deviations (10 — 15%) from scaling 
are observed in the region where the QE and A responses 
overlap. In this region effects stemming from correla¬ 
tions and 2p2h MEC can play an important role and they 
cannot be reproduced by models that assume impulsive, 
quasifree scattering on bound nucleons. Therefore these 
effects must be added to the QE and A scaling functions. 
This we do by using a parametrization na of the results 
of De Pace et al. [18| , where a fully relativistic calcula¬ 
tion of the 2p2h MEC contribution to inclusive electron 
scattering was performed. 


III. RESULTS 


QE) region, assuming that this is dominated by the A- 
resonance. We follow an updated version of the proce¬ 
dure described in (33!]. More specifically, we first define 
a non-QE experimental scaling function in this region, 
/non— qe. °pbis entails subtracting from the (e, e') double¬ 
differential cross section the SuSAv2 QE scalin g pr edic- 
tions and the 2p2h MEC contribution given by [17j : 
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Then we define a superscaling function in the region of 
the A peak as follows: 


/ non—QE (v , a ) = kp 
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where ip a is the A scaling variable and G^, G ^ are 
single-hadron functions referred to the N —> A transi¬ 
tion (see [35;] for explicit expressions). In the above om 
is the Mott cross section and vl,t are the usual kinematic 
factors. 
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In this section we first set up the model and test it 
versus electron scattering (e, e') data for the kinemat¬ 
ics relevant for the present study. Then we apply it to 
the analysis of inclusive neutrino scattering and compare 
its predictions with the data taken by the T2K and Sci- 
BooNE collaborations. 


A. The non-quasielastic scaling function 

In this subsection we construct a phenomenological 
scaling function to be used in the non-quasielastic (non- 


FIG. 1. (Color online) Averaged experimental values of 
/ non —QE(?/ j a ) together with a phenomenological fit of the non- 
QE scaling function. The colored band represents an estima¬ 
tion of the theoretical uncertainty (see text). 

Next we study the scaling behavior of J non “Q E by ana¬ 
lyzing a large set of high quality (e, e') data for 12 C, using 
similar procedures to those discussed in [l5[ • The data 
used there (see also 0 ) were chosen to match — at least 
roughly — the kinematics that are relevant for the neu¬ 
trino data under discussion; these choices of kinematics 
are listed in the figure. 
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12 C, E e = 680 MeV, 0 e = 36°, q QE * 402.5 MeV 




co(MeV) 


12 C, E e = 560 MeV, 0 e = 60°, q QE * 508.0 MeV 



12 C, E e = 560 MeV, 0 e = 145°, q QE * 795.0 MeV 



co(MeV) 


FIG. 2. Double-differential inclusive electron-carbon cross sections, da/dujdQ. The panels are labeled according to beam energy, 
scattering angle, and value of Qqe at the quasielastic peak. The results are compared with the experimental data from p~i| at 
the selected kinematics (see text). 


From this analysis, illustrated in Fig. |TJ it appears that 
scaling in the A region works reasonably well up to the 
center of the A peak, ip^ = 0, while it breaks, as ex¬ 
pected, at higher energies where other inelastic processes 
come into play. However the quality of scaling is not as 
good as in the QEP region. For this reason the non-QE 
scaling function is represented with a band, rather than 
with a function, which accounts for the spread of pseudo¬ 
data seen in Fig. [Tj This band, together with the SuSAv2 
phenomenological fits and the MEC response, can now be 
used to test the model against electron scattering data 
and to predict neutrino and antineutrino cross sections. 


B. Test versus electron scattering 

In Fig. [2] we compare the model predictions with in¬ 
clusive electron scattering data on 12 C. Although many 
high quality electron scattering data exist, here we only 
show results for a few representative choices of kinemat¬ 
ics, similar to those involved in the neutrino experiments 
that we address in the following sections. 

As observed, the model gives a good description of the 
data, provided the 2p2h MEC are included. These, as 


expected, play a major role in filling the “dip” region be¬ 
tween the QE and A peaks. The band in the final cross 
section (green region) comes from the uncertainty in the 
determination of the A superscaling function. This ex¬ 
plains that the data located in the region close to the A- 
peak are contained within the limits of the above band. 
More importantly, the model (with its associated un¬ 
certainty) is capable of reproducing successfully all data 
with particular emphasis on the dip region. This result 
gives us confidence in the reliability of the model and its 
application to the analysis of neutrino-nucleus scattering 
reactions. 


C. T2K 

In Fig. [3] we show the CC-inclusive v ^— 12 C double¬ 
differential cross section per nucleon versus the muon 
momentum, p lll for different angular bins, folded with 
the T2K flux. The QE curve corresponds to the SuSAv2 
model illustrated in Sect. HD (see [lfj for details). The 
resonant pion production curve (l7r) is obtained with the 
non-QE scaling function described above. As in the pre¬ 
vious case of electron scattering, the band corresponds to 
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FIG. 3. (Color online) The CC-inclusive T2K flux-folded i/ m - 12 C double-differential cross section per nucleon evaluated in 
the SuSAv2 model is displayed as a function of the muon momentum for different bins in the muon angle. The separate 
contributions of the QE, l7r and vector 2p2h MEC are displayed. The data are from [j. 



FIG. 4. (Color online) The CC-inclusive T2K flux-folded i/p- 
12 C double-differential cross section per nucleon is displayed 
as a function of the muon momentum, which corresponds to a 
bin in muon angle 0.00 < cos 9^ < 0.84. The full results with 
[QE+MEC+17T] and without [QE+MEC+Itt^a < 0.5)] the 
high-energy tail are shown. The data are from Tl|. 


the theoretical uncertainty associated with the extraction 
of the non-QE scaling function discussed in the previous 
section. The MEC curve corresponds to the fully rela¬ 
tivistic calculation of 2p2h excitations induced by pionic 


vector two-body current of [H} and parameterized in G2- 

We observe that the model yields excellent agreement 
with the data. Moreover, the main contribution in the 
cross section comes from the QE and pion production 
mechanisms. On the contrary, MEC play a minor role 
at these kinematics, a result that is somehow different 
from the one found in (42|. It should be noted however 
that the two calculations differ in various respects: first, 
the present model does not include the axial two-body 
current, as explained in the previous section; second, the 
two calculations, although in principle similar, involve 
different approximations in the way they account for rel¬ 
ativistic effects - the calculation of JT8( being exactly 
relativistic - and in some important technical details in 
the multidimensional integration leading to the results 
(see [H, HU ). Indeed, the MEC contributions here are so 
small that, even were AA and VA contributions that are 
as large as these VV contributions to be included, the 
net effect would still not be very significant. 

As shown in the analysis of the non-QE scaling func¬ 
tion (Fig. |T]) , scaling is not fulfilled at ip a > 0 — 0.5 due to 
other inelastic processes whose contributions start to be 
more significant at high kinematics (i.e. high momentum 
transferred). However, contributions beyond ipA = 0.5 
are not very significant at the kinematics involved in the 
v p T2K experiment, as can be seen in Fig. [fj Indeed, 
the effects in the l7r cross sections associated with this 
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p e (GeV/c) 




FIG. 5. (Color online) The CC-inclusive T2K flux-folded v e - 12 C differential cross section per nucleon evaluated in the SuSAv2 
model is displayed as a function of the electron momentum (top), cos 8 e (middle) and Q q E (bottom). The separate contributions 
of the QE, l7r and vector 2p2h MEC are displayed. The data are from Q. 


positive-'i/'A tail (high values of the transfer energy) are 
less than 10-12%. This supports the reliability of our 
model to be applied to the description of T2K muon- 
neutrino data where the average value of the neutrino 
energy is (E v ) ~ 0.85 GeV. Furthermore, upon evalu¬ 
ating the importance of high momentum transfer contri¬ 
butions by cutting the predictions at various values of 
g, we have confirmed that for muon momenta at T2K 
kinematics the model appears to be robust. 

In Fig. [5] we compare our predictions with recent 
T2K data corresponding to electron-neutrino scatter¬ 
ing 0. We show results for the flux-folded v e - 12 C dif¬ 


ferential cross section against the electron momentum 
(top panel), scattering angle (middle) and QE trans¬ 
ferred four-momentum (bottom). As observed, the model 
understimates the data at high p e and Q q E , in con¬ 
trast with the situation observed in the previous case, 
i.e., muonic neutrinos. This clearly indicates that high- 
inelasticity processes, which are not incorporated in our 
formalism, have a significant contribution in the anal¬ 
ysis of this experiment. Moreover, this is also consis¬ 
tent with the electronic neutrino flux with an average en¬ 
ergy (E Ve ) ~ 1.3 GeV, which is significantly larger than 
the value corresponding to the case of muonic neutrinos, 
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FIG. 6. (Color online) As for Fig. [5] but now only for electrons with p e > 0.55 GeV/c and cos# e > 0.72. 


(E v } ~ 0.85 GeV. This also reflects the much more im¬ 
portant tail in the electronic neutrino flux that extends 
to very high neutrino energies. From analyses of q and 
t/>A cuts as discussed above we find that effects from con¬ 
tributions lying above the A might provide 20-30% more 
strength in some cases. In particular, looking at the low¬ 
est panel in Fig. 5, one finds that the cuts have very little 
effect near the peak of the cross section at low Qq E , but 
are more significant in the tail (~30%). Thus this short¬ 
fall in the latter region is probably not unexpected. 

For completeness, we also show in Fig. [6] the “re¬ 
duced phase-space” results for the inclusive v e — 12 C re¬ 
action where a better agreement with Iow-Qq^, data is 
reached when considering only very forward processes 


(cos0 e > 0.72) and rejecting low values for the electron 
momentum (p e > 0.55 GeV/c). 


D. SciBooNE 

In this section we apply our model to the kinemat¬ 
ics involved in the SciBooNE experiment JJ, that cor¬ 
responds to CC inclusive scattering on a polystyrene 
target (CgHg). In this case the data are presented as 
a total unfolded integrated cross section as a function 
of the neutrino energy. In the unfolding procedure the 
neutrino energy is reconstructed from the kinematics of 
the outgoing lepton assuming that the process is purely 


















FIG. 7. (Color online) The CC-inclusive cross section on a 
polystyrene target (CsHs) per nucleon evaluated in the SuSA 
model as a function of the neutrino energy. The SciBooNE 
data are from [sj. 

quasielastic. This is a strong assumption, especially at 
high neutrino energies where inelastic processes become 
important, and makes the comparison of theoretical re¬ 
sults with data very delicate. We should notice that 
the averaged value of the neutrino energy for SciBooNE 
is 0.76 GeV whereas the unfolding is extended up to 5 
GeV. This makes an important difference with T2K data, 
which are not unfolded, and the corresponding calcula¬ 
tions, which involve an integration over the neutrino flux. 

Our results for inclusive scattering are compared 
with SciBooNE data in the top panel of Fig. [7] We also 
show for completeness predictions corresponding to an¬ 
tineutrinos (bottom panel). 

In the case of neutrinos (top panel), the QE contri¬ 
bution grows with the neutrino energy up to E v ~ 1 
GeV, where it saturates to a value of the cross sec¬ 
tion close to x 10~ 39 cm 2 . On the other hand, the 
resonant pion production result, l7r, (displayed as the 
brown band) starts its contribution at E v slightly be¬ 
low 0.5 GeV, and reaches its saturation value ~5 x 10” 39 
cm 2 at E v ~ 1.5 GeV. The sum of the two contributions 
agrees with the two lowest experimental points taken at 
E„ « 0.4 and 0.65 GeV, respectively. On the contrary, 
at higher neutrino energies the model clearly underpre¬ 
dicts the data by approximately a factor two at E v above 
1.5 GeV. This result clearly indicates that new channels 


and higher resonances should be included in the model 
(see also comments in the case of T2K v e data). This 
is also consistent with previous calculations based on 
RPA [42|, |44| that were restricted to neutrino energies 
below E v ~ 1.2 GeV. Although not shown in the fig¬ 
ures, we have evaluated the contribution in the resonant 
pion production ascribed to the kinematical region above 
f/>A = 0.5, i.e., the tail where high inelasticities might 
give significant effects: in this case this region (^a > 0.5) 
provides a contribution of the order of ~ 30% of the total 
integrated cross section at high energies (where satura¬ 
tion is already reached). 

Finally, as for the T2K case discussed in the previous 
section, we find the pionic 2p2h MEC to give a minor 
contribution, at most of the order of ~7% of the total 
cross section at the highest energies. The uncertainty 
introduced by the use of different parametrizations of the 
axial N —> A form factors is very small (less than 5.5% 
of the total cross section at the highest energies). 

In the case of antineutrinos (bottom panel) we observe 
that MEC contributions are about 9% at the highest en¬ 
ergies. Another important difference between neutrinos 
and antineutrinos concerns the property of saturation. 
Whereas the neutrino cross section already saturates at 
energies of the order of E V/i ~ 1.5 GeV, and this result 
applies to both QE and pionic channels, the antineutrino 
cross section continues to grow for increasing values of 
Ep^. Results in Fig . [7] are also consistent with RPA 
predictions in [42), [SJ • The analysis of T2K results on 
antineutrino CC inclusive cross sections (not completed 
yet) will undoubtedly help us in disentangling the spe¬ 
cific roles played by the different ingredients that enter 
in the description of the scattering reaction. 

IV. CONCLUSIONS 

We have compared the predictions of the recently re¬ 
vised superscaling model (SuSAv2), devised for QE scat¬ 
tering and extended to the A-resonance production re¬ 
gion, with the available inclusive data for charged cur¬ 
rent muon (electron) neutrino- 12 C reactions of the T2K 
and SciBooNE experiments, where the mean neutrino en¬ 
ergy is 0.85 GeV (1.3 GeV) and 0.8 GeV, respectively. 
The model also includes 2p2h excitations induced by vec¬ 
tor meson-exchange currents carried by the pion and has 
been tested against inclusive electron scattering. More¬ 
over, the model is fully relativistic and can therefore be 
applied at high energies, provided the relevant physics 
mechanisms are taken into account. 

Our main conclusions can be summarized as follows: 

1. The present approach provides a very good rep¬ 
resentation of the T2K experimental data. A 
similar comment applies to the SciBooNE data for 
neutrino energy below 0.7 — 0.8 GeV. 

2. On the contrary, the model fails in reproducing Sci¬ 
BooNE data at higher neutrino energies, as well as 
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T2K v e data. This is a clear signal of the rele¬ 
vance of other reaction mechanisms (not included 
in the model yet) such as resonances beyond the A, 
multi-meson production and deep inelastic scatter¬ 
ing. Work is in progress to implement these pro¬ 
cesses in the model. 

3. Pionic (vector) meson-exchange currents in neu¬ 
trino scattering are shown to play a minor role 
(<10%) for all of the kinematical situations con¬ 
sidered here. Axial-vector MEC contributions have 
yet to be included. 

4. The uncertainty related to the poorly known axial 
form factors entering in the A —>• A current does 
not present a significant impact (<6%) at the ex¬ 
perimental kinematics analyzed. 
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